Many Body Effects on the Transport Properties of Single-Molecule Devices 
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The conductance through a molecular device including electron-electron and electron-phonon 
interactions is calculated using the Numerical Renormalization Group method. At low temperatures 
and weak electron-phonon coupling the properties of the conductance can be explained in terms of 
the standard Kondo model with renormalized parameters. At large electron-phonon coupling a 
charge analog of the Kondo effect takes place that can be mapped into an anisotropic Kondo model. 
In this regime the molecule is strongly polarized by a gate voltage which leads to rectification in the 
current-voltage characteristics of the molecular junction. 

PACS numbers: 71.27.+a, 75.20.Hr, 73.63.-b, 85.65. +h 



Electronic conduction through nanoscale systems con- 
nected to external electrodes exhibits a number of special 
features that are of great importance when one considers 
their potential as electronic devices. 

The transport properties of such systems are strongly 
affected by (e-e) and electron-phonon (e-ph) interactions. 
For example, Coulomb blockade effects have been shown 
to dominate the transport in quantum dots and even 
in single molecules weakly coupled to the electrodes 
The Kondo effect has been observed in quantum dots 
and single molecules having well defined spin and charge 
states |H [H El Furthermore the effects of e-ph cou- 
pling have been observed in Inelastic Electron Tunnel- 
ing Spectra of small molecules adsorbed on surfaces Q 
and in molecular-scale transistors made of Cgo molecules 
Finally, e-ph coupling is also known to play an im- 
portant role in the transport through organic conjugated 
molecules in which transport is mediated via polaron or 
soliton propagation 

It is usually assumed that energy scales for e-e inter- 
action are much larger than those for e-ph interaction. 
However, it has been recently shown that Coulomb charg- 
ing energies of single molecules can be considerably re- 
duced by screening due to the electrodes Q. These en- 
ergies can be of the order of a few 100 meV (instead of 
a few eV for an isolated molecule) and therefore are of 
the same order of magnitude of the relaxation energies 
induced by e-ph coupling in the same systems. Interest- 
ing physics then arises when the energy scales of e-e and 
e-ph interactions merge. By studying such cases, insights 
on the origin of the features in the conductance can be 
brought and the influence of the e-ph coupling on the 
Kondo effect can be addressed. 

Although there have been many theoretical investiga- 
tions of electron transport in the presence of either e-e or 
e-ph interactions, the case in which both interactions are 
present has only been studied recently 0, 0, either 
in the high temperature or weak e-ph coupling regimes. 

In this Letter, we report results of non-pertubative 
calculations of the linear transport properties through a 
molecular device including e-e and e-ph interactions. The 
calculations are performed using the Numerical Renor- 



malization Group (NRG) technique for a broad range of 
characteristic parameters. 

We study a model of a molecule with a single relevant 
electronic level coupled linearly to a vibrational mode of 
frequency ujo and to the left (L) and right (R) metallic 
electrodes. The Hamiltonian of the system is 



H = H 



M 



H E + H 



M-E, 



(1) 



where the first two terms describe the isolated molecule 
and electrodes, respectively, and the last term describes 
their coupling. We have 

H M = e d n d + Un d ^n dl - A {n d - 1) (a + a f ) (2) 



-ojQa' a 

E 



e a (k) c 



kcra kcra ' 



k,er,a={L,R} 
H M -E = Vka { d ° C k™ + C icra d cr 



(3) 

(4) 



Here, n d = ^t^a is the charge of the molecule, e d is 
the position of the electronic molecular level relative to 
the Fermi level of the electrodes and U is the Coulomb 
repulsion between two electrons that occupy the same 
molecular level, creates a phonon of frequency ojq 
and A is the e-ph coupling constant. We consider for 
simplicity the case of identical electrodes and contacts. In 
the wide-band limit the conductance G of the molecular 
junction at zero bias is (l^fl4| 



dl 



2e2 a 

= 7rA 
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duj — 



dm 



Pdiuj) , 
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where p d (w) = —tt 1 ImG dd (uj), G dd (ui) is the exact 
electronic Green function of the molecule in the pres- 
ence of the leads, f(oj) the Fermi distribution and A = 
2ttpo(V£)fs- Here, p is the electrodes' electronic density 
of states at the Fermi level and the brackets denote an av- 
erage over the Fermi surface. Hamiltonian is electron- 
hole symmetric for all values of A at e d = e* d = — U/2 
where (n d ) = 1. This point we will referred to as the 
symmetric point in the following. 
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At zero temperature Eq. JSJ) reads 

G = G TrA Pd (0) = Go sin 2 {^{n d )/2) , (6) 

where Go = 2e 2 /h is the quantum of conductance and the 
second equality follows from Luttinger's theorem ^3- At 
T = the conductance G thus takes its maximum value 
Go at the symmetric point. 

Before presenting the results of a full numerical solu- 
tion of the problem we discuss two limiting cases that 
can be treated analytically providing us with a qualita- 
tive picture of the dependence of G upon T and £ d - 

The Hamiltonian Hm of the isolated molecule can be 
readily diagonalized. The eigenfunctions are direct prod- 
ucts of electronic states (denoted by a subscript e) and 
oscillator states. The eigenstates and eigenvalues are 



\Q,m) = U-\0) e \m), 
\a,m) = \a) e \m), 

\2,m) = U+\n) e \m), 
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+ 2s d + U + mu> , 



where U ± — exp [±A/cjo (a* — a)] and \m) is the m-th 
excited state of the harmonic oscillator. Two limiting 
cases can be considered. 

Weak electron-phonon coupling, 2\ 2 /ujq <C U. In this 
case the ground-state of the isolated molecule is the spin- 
doublet | (7,0). There is a large charge excitation gap 
~ C/ e ff = U — 2X 2 /loq at the symmetric point and the 
low-energy excitations of the full system are spin fluctua- 
tions described by the usual Kondo Hamiltonian |l5j . Us- 
ing standard second-order perturbation theory the cou- 
pling constant can be expressed in terms of a single 
matrix element: (g,Q\Hm— E\0,'m)(0,m\HM-E\o',0) = 
(a,0\H M -E\0,m){0,m\H M -E\ -ff,0) cx \{0\U+\m)\ 2 . 
At Ed — £d we find the Kondo coupling constant [la] 
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where |(0|J7 + |to)| 5 



-( A /"») 2 (A/w ) 2m /m!. 



Below the Kondo temperature Tk oc exp [—\/{JkPo)], 
G(£* d ) ~ Go- Expanding Eq. JSJ around A = we find 
Jk/Jk(0) pa 1 + 2 [1 + U/{2u )r 1 \ 2 /{Uuj q ). The e-ph 
coupling thus leads to an increase of the Kondo temper- 
ature in this regime. For T > Tk, Coulomb blockade 
peaks separated by U c s are expected at gate voltages 
E d PS -A 2 /w and s d ~ -U + \ 2 /u> . 

Expanding Eq. (jSJ) around the symmetric point we 

have G(£ d ) - G{e* d ) G (x c 7r/2) 2 (e d - £* d ) 2 , where 

Xc is the charge susceptibility, inversely proportional to 
the charge gap. The width of the T = conductance 
peak, Ae^ pa Xc 1 ~ U c g, decreases with increasing A in 
the weak coupling regime. 

Strong electron-phonon coupling, 2X 2 /ujq U. In this 
regime the ground-state doublet of the isolated molecule 



is composed of the states |2,0) and |0, 0), degener- 
ate at the symmetric point. The low-energy excita- 
tions of the full system are now charge fluctuations and 
there is a large gap for spin fluctuations. The low- 
energy excitations are described by an effective Kondo 
model in which the role of the spin is played by a 
pseudo-spin variable that represents the two states of 
the lowest lying doublet. A priori there is no ro- 
tational invariance in pseudo-spin space and, in fact, 
in the effective Hamiltonian two different matrix ele- 
ments appear: (2, 0\Hm-e\<7, Tn)(a, m\HM-E |2, 0} oc 
|(0|^+|m)| 2 , and {2,Q\H M _ E \a,m){a,m\H M ^ E \Q,0) cx 
(0\U-\m)(m\U + \0) = (-l) m |(0|C/ + |m)| 2 . The effective 
model is thus the anisotropic Kondo model (AKM) with 
couplings Jii and J± given by 
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and Kondo temperature Takm given by |l7j : 

Takm <xD(j±/J n )^ . (10) 

Asymptotically, J±/Ju ~ exp [— 2(A/w ) 2 ] and Takm ~ 
exp [— a(A/wo) 4 ] with a — ttujq/A. In this regime the 
Kondo temperature decreases sharply with increasing A. 
In contrast to the weak e-ph coupling case, no Coulomb 
blockade peaks are expected for T > Takm- Another 
important difference from the weak A case is that, now, 
the charge susceptibility (cx ?akm) ^ s ver y large. Then, 
the width of the conductance peak versus gate voltage 
As d pa Takm also decreases sharply as A increases. The 
region in the T — Ed plane where G pa Go is thus strongly 
suppressed by a strong e-ph coupling. 

We now turn to the discussion of our numerical re- 
sults obtained using the NRG method [3 EH]. The 
original NRG method has been modified to include e-ph 
coupling |2fl | and to calculate the spectral density Pd(to) 
according to Ref. [U- We have performed calculations 
for a broad range of parameters. We set k E = 1 and 
take the half-bandwidth of the electrodes as the unit of 
energy. In the following, we present results for the set of 
parameters A = 0.016, w = 0.05 and U = 0.1. 

Figure n shows the linear conductance G as a function 
of the gate voltage Ed for T = and T — A and several 
values of the e-ph coupling A. For A = jFig Uall and 
T = A the Coulomb blockade peaks [13, separated 
by the charging energy U are clearly seen. With decreas- 
ing T the conductance in the Coulomb blockade valley 
increases as the Kondo effect develops. At zero temper- 
ature G has a single peak of width Asd ~ U centered at 
E* d . The main effect of a weak coupling to the vibrational 
mode at finite T [FigOJb)] is a reduction of the distance 
between the Coulomb blockade peaks, now given by C/eff 
as anticipated above. The width of the T = peak is also 
of the order of U e s . This large width is a consequence of 
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FIG. 1: Conductance G versus gate voltage Ed for different 
values of A (a-d). (thick lines: T = 0, thin lines: T = 0.016). 
Parameters are U = 0.1, A = 0.016 and uio = 0.05. 
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FIG. 3: Full width at half maximum F of the central peak 
in the spectral density pd versus A at the symmetric point 
Ed = £*d- Solid lines are fits based on the given analytical 
expressions (see text). Inset: Pd(u) for three values of A. 
Other parameters as in Fig. 
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FIG. 2: Zero-bias conductance G versus temperature T for 
different values of A. (a) Symmetric point Ed = —U/2 — 
—0.05. (b) Asymmetric case, 5£d = Ed — E* d = 0.01. Other 
parameters as in Fig. Q 



the rigidity of the ground state of the molecule against 
charge fluctuations. 

For strong e-ph coupling [Fig^d)] the features change 
qualitatively. A single peak is observed in the conduc- 
tance at all temperatures and its width sharply decreases 
with T . The peak narrowing results from the dramatic 
increase of the charge susceptibility of the molecule in 
the strong coupling regime. Moving away from the sym- 
metric point by application of a gate voltage produces 
an effect similar to that of applying a magnetic field in 
the standard Kondo effect |2J|: the degeneracy of the 
ground-state doublet is broken and the Kondo resonance 
in the spectral density is destroyed |2fij . 

For U e g < A [Fig^c)] the system is in a mixed valence 
regime in which the four charge states of the molecule 
are nearly degenerate for Ed ~ £* d - A single peak is then 
observed in the conductance at all temperatures with a 
T = width that is determined by the hybridization A. 

In all regimes discussed so far, the T = conductance 
at e* d is perfect as required by Eq. |JBJ. 

Figure |21 shows the T-dependence of G for several val- 



ues of A. At the symmetric point [Figl^a)] G ~ Go 
below a temperature scale that depends on the coupling. 
This scale varies non-monotonically with A, first increas- 
ing and then decreasing. For very small A the behavior in 
the asymmetric case [Figffib)] is similar to the previous 
one. For the larger couplings, however, G < Go at all 
temperatures. 

The inset in Fig. [31 shows the spectral density Pd{u) 
at the symmetric point for three values of the e-ph cou- 
pling. The width of the central peak is renormalized by 
the interaction. It is also non-monotonic as a function of 
A as shown in the main plot where fits to the theoretical 
expressions discussed above are also shown. For weak 
coupling region we use T = A exp [—\/{JkPo)\ with Jk 
given by Eq. © . For strong coupling we use Y — A' Takm 
with Takm given by Ea. HlUI) . instead. Only the ampli- 
tudes A, A 1 are fitting parameters. It can be seen from 
Fig-OHthat the agreement with the predicted widths F(A) 
is excellent. For intermediate couplings the charge fluc- 
tuations are no longer blocked by the Coulomb interac- 
tion and the system is in the mixed valence regime where 

r = 0(A). 

Figure 0] shows the spectral density for strong e-ph 
coupling and several values of 5ed < 0. At the sym- 
metric point there are peaks at the energies of the iso- 
lated molecule given by Eq.Q. There is also a peak 
at the Fermi level associated with the charge Kondo ef- 
fect. Its width is proportional to Takm [Eq. l[TU)l]. For 
1 8sd | -C Takm there is little change in the spectral den- 
sity. When the asymmetry increases, spectral weight is 
transferred from positive to negative energies ui and Pd(0) 
is suppressed. For \5ed\ > Takm, the central peak disap- 
pears and almost all the spectral weight is concentrated 
in the negative oj region. The behavior for positive Sed 
is obtained through the transformation lu — > — u>. This 
spectral polarization is the result of the strong charge 
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junction are controlled by the large charge polarizabil- 
ity of the molecule. The latter induces a strong spectral 
polarization that leads to a rectifying behaviour of the 
current-voltage characteristics of the molecular junction. 

We thank C. Castellani, M. Grilli and M. J. Rozenberg 
for useful discussions. 



FIG. 4: Spectral density pd for different values of the asym- 
metry 5ed = Ed — Ed in the strong e-ph coupling regime 
(A = 1.2a>o). Here Takm ~ 0.0012. Left inset : average 
charge (rid) in the molecule as a function of the asymmetry. 
Right inset: current- voltage curves showing rectification in 
the asymmetric case. Other parameters as in Fig. Q 



polarization induced by the asymmetry. It is illustrated 
in the left inset to Figure 0] that shows the ground-state 
charge of the molecule (n d ) as a function of 5e d . It can be 
seen that charge fluctuations in the molecule are rapidly 
suppressed for \6ed\ > Takm- 

An interesting consequence of this spectral asymme- 
try is that the molecular device exhibits strong rectify- 
ing properties in the strong e-ph coupling limit. The 
right inset to Fig. 21 shows the current-voltage (/ — V) 
characteristics for the device in this regime |26|. At the 
symmetric point, the I — V curve is symmetric around 
zero bias V and it exhibits some structures associated to 
the peaks in the spectral density. In the asymmetric case 
(with Sed < 0), the molecule is polarized and the I — V 
curve is asymmetric: the current for negative bias V is 
much larger than for positive V (the opposite occurs for 
Se d > 0). 

In summary, we have studied the transport properties 
of a molecular device including e-e and e-ph interactions. 
We have shown that at low temperatures and weak e- 
ph coupling, the conductance properties of the molecu- 
lar junction are described by the standard Kondo model 
with renormalized parameters. At large e-ph coupling, 
new physics appears: the properties of the device are de- 
scribed in terms of a charge analog of the Kondo effect. 
In this regime the transport properties of the molecular 
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